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1.0 Quality Compression Optimization 

.wi°I 10 n Bg 'a 110,18 ^ With ^ ira P rovem «ot of the compressed 
mage quality As stated m the report introduction, we believe th" is 
the last significant major problem to be addressed in the still 
imagery compression system. 



1.1 Introduction 



The ' issue is the mechanism which can be used to improve the quality 
of the im age based upon the nature of the parameters used in the 
compression process. We think this is an important issue, and in our 
expenmemation. we have found fairly striking results. For example 
if we hold the compression ratio constant, we can compress an image 
and obtain on a subjective scale of C. where 

• A is no observable defect, 

• B is observable but not noticeable, 

• C is quite noticeable, but not distracting from the image 

• D is very noticeable and detracts from the image, 

• E is unacceptable. 

If we now change the parameters by a hand optimization, we find we 
can change the subjective evaluation from a C to a B with the 
compression ratio left alone. We believe this is significant, in that it 
guarantees a compression system which is tailored for each image 
independently, rather than have each image compressed by the same 
set of parameters irrespective of the image content. 

Figie Vi 49 Pr ° CCSS ° f COmprCSsion ' consider following figure. 
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Figure 49 

RGB Image in 24 Bit Color Depth Showing Transform to YIQ 

Figure 49 shows the relationship between the RGB and YIQ color 
representation and the processes which take place in the algorithm 
as this data is transformed into the final lossy set of wavelet 
coefficients.. The processes which apply to the imagery are shown as 
heavy arrows with the process identifier shown as part of the arrow. 
If no process number is present, then the arrow is just a passive link 
between processes and data. The important point to note is that in 
color imagery, we deal with three images (one luminescence, and two 
color planes). The Y component is critical, while the IQ components 
are less sensitive to error introduced by the compression system. 

The processes shown in Figure 49 are as follows: 

• (1) Transform the RGB format into YIQ format, using long 
integers. This format is only an approximation of the YIQ format. 

• (2) Transform the YIQ planes into a wavelet decomposition using 
our own integer wavelet transform. This produces the result 
shown in Figure 50, but for each plane. 

• (3) We have shown an alternative process which is an optional 
down sampling for the IQ color planes. This down sampling may 
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be done once or twice to produce two image planes either one 
fourth or one sixteenth the size of the original plane. If this 
process is to be done, it will be accomplished prior to the wavelet 
transform of process (2). 

• (4) Here the first step in the loss occurs. The wavelet coefficients 
are now quantized to a number of levels depending upon which 
quadrant is being processed, and the desired compression or 
quality factor. 

• (5) Simultaneously with step (4), the wavelet coefficients are 
being matched against threshold values, and if the values are less 
than the established threshold values specified, then the resultant 
value is set to zero. 

• (6) The last step in the process is to entropy compress the 
resultant coefficients using cither Arithmetic, Run Length, or 
Huffman, of Huffman and Run Length combined. The key issue is 
the amount of compression desired against the invested time 
required to get that level of compression. 

The issue now, irrespective of the down sampling or not of the IQ 
components, is the fact that we process the three planes into five 
levels in a decomposition as shown in Figure 50. From this figure, 
one can see a total of 16 regions (quadrants) which are defined by 
the numbers 1 through 16. Each of these sixteen quadrants have two 
additional parameters associated with them. The parameters define 
the quantization and threshold values for that particular quadrant. 
Since there are three planes for color (only one for gray level) the 
maximum number of parameters that the user can control is 96 -- 16 
for quantization and 16 for thresholding for each of the three color 
layers. In the case of a gray level image, there are only parameters 
for one layer. 

Our experience has shown the parameters for an image are very 
sensitive in some cases, and not in other cases. In order to measure 
this sensitivity, we generated the variance of the wavelet coefficients 
in the 16 quadrants. Table 1 provides these values for each of the 
three planes. 
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Figure 50 

Decomposition of a Plane into Five Levels 
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56 
23 
10 
128 
55 
37 
253 
75 
64 
499 
138 
156 



Tulips 
Wavelet Varia^ K 



4.7 
3.2 
1.3 
5.8 
4.5 
2.8 
11.4 
7.4 
3.3 
27.5 
11.8 
6.5 
41.7 
11.2 
12.2 
1-281 8 



1.4 
1.1 
0.6 
1.7 
1.1 
0.8 
2.6 
1.1 
1.0 
4.1 
1.5 
1.1 
10.1 
2.5 
2.0 



164.2 
134.6 
28.9 
276.1 
274.3 
124.8 
225.7 
298.1 
114.2 
174.8 
285.1 
107.9 
135.4 
245.0 
89.3 
7-408* 



_L 



27.3 
24.9 
4.6 
57.7 
67.3 
22.7 
73.2 

107.3 
35.5 
73.9 

128.2 
47.8 
75.6 

144.0 
49.1 

690.4 



6.6 
5.0 
1.1 
12.9 
12.4 
5.0 
13.7 
16.2 
6.8 
10.7 
14.9 
6.7 
8.7 
11.2 
5.1 
77.9 



Table 1 

Wavelet Coefficient Variance by Numbered Quadrant for Two Images 
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The clear information to be gained by the numbers in Table 1 is the 
images are quite different, and the quantization or threshold levels 
set for all images will only be an approximate solution at best The 
optimum solution would be to have the quantization and threshold 
values set according to the variance values. Such settings could be 
found in a table with various ranges, and for each such range, the 
parameters of interest could be defined. This would give a more 
optimal solution, but still not the optimal solution. In order to get 
optimality, one would need to search over the variable space using 
the near optimal settings to find the actual best values for the 
parameters. 

As the values in Table 1 grow, the implication is a finer mesh size for 
quantization, That is, the number of bits one wishes to allocate to the 
output could be varied by quadrant. Those quadrants with large 
variances will utilize more bits, while those with low variants will 
utilize fewer bits. In this way, the number of bits resulting from 
quantization will remain the same, but their allocation will differ 
depending upon the nature of the image. 

1.2 Approach 

The solution to the problem posed in the previous section is to 
determine how the ninety-six parameters interact with one another. 
The problem is a bit more sophisticated than just measuring 
parameters, however. The issue is with each parameter change, the 
compression ratio will change. If a compression ratio, or a quality 
factor which indirectly defines a compression ratio, is specified, then 
the user wants the compression ratio to remain identical over the 
changes in the parameters. In order to accomplish this, there are 
two parameters which we must monitor: PSNR (peak signal to noise 
ratio which is defined to be PSNR = 20 logio (X/MSE) where the X is 
the average absolute value of the pixels in the after image and MSE 
is the mean squared error measured between the before and after 
image) and the compression ratio. The compression ratio must be 
held constant, and the PSNR needs to increase, and the way to 
increase the PSNR is to reduce the MSE, 

The difficulty with this system as described is in many cases, small 
changes in the parameters introduce significant changes in the MSE. 
Also, we believe, the parameters are not independent. We have also 
seen images where the parameters can be changed in one way, then 
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altered, and the results are exactly the same. This indicates the 
optimal value is not a single point, but rather something like a plane 
with little slope. 



1.3 Status 

We have established the rules under which the optimal solution will 
need to exist, and arc at the moment writing software to measure the 
variance within the Lightning Strike environment. Once this is done, 
we will being examining many images to sec how close we can 
determine the optimal parameters for a defined set of variances. 
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CIS-2 Image Compression Algorithm 



HONGYANG CHAO 



Part I: Brief Review of LSIC 3.0 



1. Introduction 

CIS-2 (temporary name), which has been being used in Lightning Strike 3.0 image 
compression software, is a wavelet based image compression algorithm. ClS-2 has 
following inventions: 

• Integer reversible wavelet algorithm with Property of Precision Preservation; 

• Subband oriented quantization and related entropy coding; 

• Wavelet lossless compression for color and gray images; 

• Progressive transmission algorithm for color bit compression; 

• Progressive transmission and decompression algorithms; 
9 Non-unifonn image compression algorithm; 

• Quality based wavelet coefficient quantization tables; 

• Attached optional post-processing filters; 

• Image map editor; 

• Optional peak signal noise ratio controlled compression; 

• Special split and merge wavelet compression algorithm for very big image 
compression without any boundary effects ; 

• Image dependent parameter optimization • 

2. Main steps of the algorithm 

In Lsic 3,0, three kinds of different image compression methods arc included: 



Section 2*1-23 will give brief description of above method The details will discuss 
later. 

2.1, Main step of method 1 

Figure 1 and Figure 2 give the flow charts of the image compression and 
decompression of methodl respectively. Every step in both compression and 
decompression has lot of details, which will be described later. 



Method 1 : Quality controlled wavelet based compression 
Method 2: Color bit depth compression 
Method 3: Wavelet lossless compression 
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Figure 1: Compression flow chart for Method 1 
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Figure 2: Decompression flow chart for Method 1 
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22. Main step of method 2 

This method based on using less number of colors to approximately represent 
original images. Following figure gives the main step of the algorithm for the 
compression and decompression. 



Input Original Image 
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Figure 3: Compression flow chart for Method 2 
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Figure 4: Decompression flow chart for Method 2 

23, Main ateps of method 3 

The main steps of method 3 is almost as same as method 1. However, at the every 
step we use different methods. Following are its compression and decompression 
flow charts: 
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Figure 5: Compression flow chart for Method 3 
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Figure 6: Decompression flow chart for Method 3 



3* Brief Description for Other Features 

Most Features (or inventions) are included in above three compression methods. 
Following is the brief introduction for some inventions list above. 
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This algorithm allows users to decode images from the lowest resolution to 
highest resolution. The advantage of this feature is that users can download small 
piece of the coded file and view the image at lower resolution to determine if they 
want to download the whole image. 

Steps: 

(a) Input the lowest pass component LL° of the coded file and reconstruct the 
lowest resolution image 7° ; 

(b) Display image 7° . If the user doesn't like it or the resolution is big enough for 
, stop; otherwise, go to next step; 

(c) Input the lowest three band-pass components HL\ LH'and HH° 
successively in the current coded file. Reconstruct the new image 7' from 
LL\ HL\ LH° and 7/77° . Let 7° -7 1 , go to step (b). 

3.2. Non-uniform image compression 

The algorithm allows users to divide the image into several parts with different 
interests and compress these areas with different qualities. The areas can have any 
shape. This algorithm is only available for method 1, 

Original image goes though all of the procedure except quantization part, which 
follows the steps below; 

(a) Creating the bitmap matrices related to the areas chosen by the user; 

(b) Wavelet transform to every bitmap matrix; 

(c) Different quantization in different areas according to the transformed matrices 
obtained above step. 

33. Peak Signal Noise Ratio (PSNR) controlled compression 

Peak Signal Noise Ratio (PSNR) is an image quality measurement used by most 
professional people. PSNR controlled compression allows users to choose their 
desired PSNR for the compressed image. 

The related algorithm is an iterated system: 

(a) Picking an initial parameter setting ? 0 ; 

(b) Quantize the wavelet coefficients with /> and calculate the corresponding 



(c) If the PSNR is close to desired one, stop and output the coded file; otherwise, 
get an adjusted vector A/> and set P 0 <- /> + AP Q , go to step (b); 



PSNR; 
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3.4. Attached optional post-processing filten 

Users can choose any number of following processing filters at their compressing 
time. The desired results can stored in the coded image file, and, anyone who 
decompress the coded file will see the same result immediately. 

• Sharpening images 

• Smoothing images 

• Improving the visual quality 

• Brightening the images 

3*5. Image map editor 

Image Map Editor create a image map over Lsic3.0 compressed image file. User 
selects one or several areas of compressed image, assigns the http links to the areas, 
then. Image Map Editor calculates the coordinates of the areas and outputs a HTML 
associate with the image. User can add such information into his/her source code. 

Fol lowing is an exa mple of such image map: 
ycflSdBED SRC-^co^reoii^type = "image/cis-cod" WIDTH* "257" poly= "44, 
^2> V 4S, 103, 78, 103, 86, 54, 86, 54^^href = "http://www.infinop.com: , X^MBED> 

3.6. Split and merge wavelet algorithm for very Wg image compression 

This algorithm allows users to compress very big image by an ordinary machine. 
The key is to divide the original image into several smaller pieces and 
compress/decompress them separately by using overlap and disoverlap technique. 
With this technique, the compression/decompression piece by piece is equal to 
compress/decompress the whole image together, which means users won't see any 
edge effect at decompressed image which appears at general split method. 

Also, with this algorithm, users can either decompress the whole image or choose 
the specific part to decompress according to a image map we create for the division . 

3.7, Image dependent optimized parameter setting 

This algorithm allows user to get the best (or almost best ) image quality at the 
desired compression ratio by choosing image related parameter setting. 

3 A Integer reversible wavelet algorithm with PPP property 

See attached unpublished paper titled as 41 An Approach to Fast Integer Reversible 
Wavelet Transformations for Image Compression" 
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3.9. Integer color transformation 

This algorithm is a integer reversible transform which has been used in lossless 
color image compression (Method 3) for Lsic 3.0. 

The algorithm transform RGB color components to a new set of color 
components Y-Nb-Nr: 

Forward transform RGB to Y-Nb-Nr: 

r = G + Int(Y). 
M>«J-Int(f), 

Afr = tf-Int(f). 
Inverse transform Y-Nb-Nr to RGB: 

*-*r + lnt(|), 
S = A& + Int(f), 
G=r-Int(^). 

3.10* Sabband related Quantization and entropy coding 

This entropy coding method is just designed for wavelet based image 
compression. The main idea is to take the advantage of different quantization at 
different subbands and encode each band according to its content This method 
reduce the coding cost greatly. 



CIS-1 Image Compression Algorithm 

HONGYANG CHAO* 
Computer and Information Science Inc. 

1. Introduction 

CIS-1, which has been being used in Lightning Strike image compression software, is a 
wavelet based image compression algorithm. CIS-1 has following advantages: 
o Reach almost optimal compression ratio; 

o Keep the major characteristics as more as possible. In other words, it reduce 

insignificant components gradually according to human visual system, so that people 
can still accept the image quality at the extremely high compression ratio; 

o Fast 

2. Main steps of the algorithm 

Figure 1 and Figure 2 give the flow charts of the image compression and decompression 
respectively. Every step in both compression and decompression has lot of details, which 
will be described later. 

Input Original Image 
Display Image 
Image Tran sformations 

r 

Quantization an d Thresholding 

En coding 

Output Coded Image File 
Figure 1: Image compression 
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Figure 2: Image decompression 

3. Image Compression 

For the compression, we only have to describe three parts: image transformations, 
quantization and thresholding, and entropy coding. 

3.1. Image Transformations 

The image transformations involved in this algorithm include color transform (for color 
images) and wavelet decomposition ( for both gray level images and color images). 

(a) Color transform 

In general, input color images are based on RGB color model, such as 

TIFF or BMP images. In order to get high compression ratio, it is better to change RGB 

color model to other color models, such as YIQ or YUV models. 

RGB to YIQ: 

IYI R>.299 0.587 0.114 T*l 
1/ l = lo.596 -0.275 -0.321 1 G I. 
LgJ Lo.212 -0.523 0.311 JjfJ 

RGB to YUV: 



["0.299 0.587 0.114T*! 
I U\ = 1-0.148 -0.289 0.439 I G I. 
[v\ Lo.615 -0.515 -0.1 1*J 



(b) Wavelet decomposition 

The purpose of wavelet transform is to represent the original image by different basis to 
achieve the objective of decorrelation. There are a lot of wavelets which can be used in 
this step. In CIS-1, we use a wavelet which results in the following algorithm: Suppose 
C ° = [ c *l** (/-O/'^A/" 1 ; *-0,'",tf-l) is original image, where M and// are 
integers which have the common factor 2 L (L is a positive integer). After one-level 
wavelet decomposition, we will get four parts as shown in figure 3. 
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Figure 3. Wavelet image decomposition 

We call C 1 - [c^] (; = 0, • • ¥ - 1 ; k = 0,- - y£ - 1 ) the blurred image of C° , HD 1 the 
horizontal high frequency part of C°, VD 1 the vertical high frequency part, and DD 1 the 
diagonal ones. Setting C° = C 1 , we can repeat the same procedure L times or until the 
size of the new blurred image C 1 is small enough. Therefore, we only have to give the 
algorithm for one level decomposition: 

(1) Let C° = rC° , r >0 is a factor which can be changed for different needs. 

(2) Transform for image columns: 
o For ifcsO/ calculate 



*0 *0 



o For Jt = 0,- -,N-l, calculate 



(3.1.1) 



-f \ 



} 




(3.1.2) 



(3) Transform for rows: 

o For y :=(),••.,¥ -1, computing 

>«/i t« 



M J.o = — — . 



and 



(3.1.3) 



|*vo ~ c j,\ 2 » 

yjk=c jaM '* 2 * = l,...,¥-2, (3.1.4) 



o For ; = 0, • ^ - 1 , computing 



and 



(3.1.5) 



| va /> -<•/.! j • 

jv4-^,-^^,it = l,.^-2, (3.1.6) 

(4) C = [c^], HD' = VD' = [v^Jand DD 1 = [d^} y = 0,-,f -1 

* = 0,-,f-l. 

Remark. If it is necessary, we also can use matrix multiply 

Wavalet Coefficient Image of / levels = W/C°Wj. 



9>->l> 



Here, W, is the transform matrix for / level wavelet decomposition. 
3.2. Thresholding and Quantization 

Both thresholding and Quantization allow us to reduce accuracy with which the wavelet 
coefficients are represented when converting the wavelet decomposition to an integer 
representation. This can be very important in image compression, as it tends to make 
many coefficients zeros-especially those for high spatial frequencies. 

After L level, for example 1=3, wavelet decomposition, we get the wavelet coefficients 
of the original image as plotted in Figure 4: 
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Figure 4. L»3 wavelet coefficients distribution 

(a) Thresholding 

In algorithm CIS-1, we use multilevel uniform thresholding method: Let 
u u , T -«.-.'*.'*i> 

be the chosen thresholds, where/, is the threshold for / th (/ = 1, • • -,£,) level and / £+1 is a 
threshold for blurred image . Thresholding is to set every entry in the blocks &, 



HD' , VD' and DD' (/ = 1, • . ,L) to be zeros if its absolute value is not greater than the 
corresponding threshold. 

Remark For color image, we can have three threshold vectors which correspond three 
different color bands, such as Y, I and Q. 

(b) Quantization 

Quantization is to scale the wavelet coefficients and truncate them to integer values. In 
CIS-1, we use the quantization table shown in Table 1 to implement it 





<Ikd 
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Table 1. Quantization table 
Here, the entries q ^ are quantization factors for blocks HD' (/ = 1, • • , and q 
for blocks VD' and DD ; (/ « 1, . . -,Z) respectively, and the factor q? 1 is for the most 
blurred image C 1 . All the factors are integers between 0 and 255. The quantization 
scheme for the block HD' (/ - 1, ••-,£) is 

^ =ro ^(-^r).;=l>."^-l; *«0,-».#-l. (3.2.1) 



I 

DD 



max. 



Here, hcfj* (j = 0, — , $ - 1; k = 0,- • ;f - 1 ) are quantized wavelet coefficients in block 
HD ; (/ = 1,-,I), 



max™ = max 

and the function round (x) gives the nearest integer of x . The scheme of quantization for 
other blocks are the similar to (3.2.1). 

Remark For color image, as the same as thresholding, we can have three separate 
quantization tables for different color bands. 

3 J. Entropy Coding 

Here, the encoding means the lossless compression for the wavelet coefficients. It is 
divided into two parts: Coefficient arrangement and entropy coding (Huffman or 
arithmetic). 



4. Decompression 

4.1 Decoding 

Decoding, just as encoding, can be divided into two parts: Entropy decoding (Hufiman or 
arithmetic), and coefficient rearranging. 

4.2 Dequantizatloa 

After Decoding, we get quantized wavelet coefficients in 3 *L+7 Blocks, Dequantizing 
uses the same quantization table as quantizing, and the scheme as follow: for / = 1, • • \L 

<Ihd 

(4.2. 1) allows us to get the approximate coefficients for the blocks HD' (/ = !,- • , 
which is shown in Figure 4 . The dequantizing scheme for other blocks are similar to 

(4.1.2) . 



43 Inverse Image Transformations 
(a) Wavelet reconstruction 

We are going to describe the algorithm for one-level reconstruction which is plotted in 
Figure 5. 
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Figure 5. One-level wavelet reconstruction 

(1) Inverse transform for rows: 
o For /sO,"-,^-!, calculate 

+ k = \,..A-2, (4.3.1) 

I 

and 



o For j = O,"-,^ - 1 , calculate 

i 2 >.i !3< vo + — 2 — 



y - 



^, *-l.-.4f-2. 



1 =i 



\ C J.»- 



= c 



■hct. 



(4.3.2) 



(4.3.3) 



and 



jo fl ina j.o' 
1 



(2) Inverse transform for column: 

o For k = 0,- • N - 1 , calculate 

f 



and 



2 ' 



I 



- 0 _ -1 JJI 



(4.3.4) 



(4.3.5) 



°- 20 -2* 



<1* 



'0*» 



c° = - (c +c° ^-2J' / = 1 L 

-ft''. y=o,-,A/-i; * = o,-,^-i.c o =[c; t ^. 



(4.3.6) 



(b) Inverse color transform 

For color image, we have to do inverse color transform 
o YIQ to RGB 

[R] ["1 000 0.956 0.621 JY~[ 
\G 1 = 1 1.000 -0.272 -0.64711 /I. 
UJ Ll-000 -1.106 1.703 JLfiJ 



'X. 



(4.3.2) 



o For y = 0,- ••,^-1, calculate 



2 

Ad! i + Adli 



•JM+\ ~ c Jk + 0 t * - i» ,T 



(4.3.3) 



and 



JO C J\ * nU J,Q> 



1 



(2) Inverse transform for column: 

o For k = 0,- calculate 



+ c! M . 1 )-2W' * = l,-f-l. 



(4.3.4) 



and 



(4.3.5) 



[ 



«o°*"^'~ 2 4*» (4.3.6) 

(3)c; jk =c; <t /r, /=0, -,M-1; *-0,.-.JV-l.C , -[cJ t ] fii1 . 
(b) Inverse color transform 

For color image, we have to do inverse color transform 
o YIQ to RGB 

1**1 ["1.000 0.956 0.621 1X1 
I G 1 = 1 1.000 -0.272 -0.64711 / I. 
[b\ b-000 -1.106 1.703 JLfiJ 



o YUVtoRGB 

UJ 



4.4 Necessary Pott Image Processing 

(a) Color Quatization 



ri-000 0.000 L140TI1 
1 1.000 -0.395 -0.581 If/ 1 
Ll.000 2.032 0.000 J. kJ 
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Abstract 

In this paper, we propose a general method for creating integer wavelet transforms which can be 
used in both lossless (reversible) and lossy compression of signals and images with arbitrary size. 
This method allows us to get a series of transformations which are very close to the conesoonding 
biorthogonal wavelet transforms or some non-orthogonal wavelet transforms, but can be calculated 
with only integer addition and bit-shift operations. In addition, the integer wavelet transforms 
created in this paper possess a property of precision preservation (PPP). This property is very 
useful for conserving memory in both compression and decompression, and speed up the whole 
procedure in some applications. The motivation of this paper comes from the lifting scheme [1] 
and S+P transform [3j. 
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1. Introduction 

The wavelet transform has proven to be one of the most powerful cools in the field of image 
compression. Theoretically, the wavelet transformation is lossless, but since all computers have 
only finite precision, most of transformations are lossy in practice, even when we use floating 
point calculations. On the other hand, integer calculations are much faster than floating point for 
virtually all computers; and integer computations are much easier to implement in hardware which 
is more important in some applications. The memory utilization of integers is also a positive 
consideration. The difficulty is, if we directly use integers in the wavelet transform and its inverse 
without some proper considerations, it will cause the loss of accuracy. For some important image 
applications, the user wants to have complete control of the precision in which the image pixels are 
represented during the compression process, and thus prefers to have the image compressed from 
lossless to lossy. 

Lossless compression is also very important for images found in such applications as 
medical and space science. In such situations, the designer of the compression algorithm must be 
very careful to avoid discarding any information that may be required or even useful at some later 
point. From the academic point of view, it is also very interesting to have a compression scheme 
which has very fast performance, and which can exactly reconstruct the image when necessary. In 
addition, it is also very useful to have some wavelet transforms which exhibit the property of 
precision preservation (PPP), which can be utilised in the computer which has limited precision 
and limited memory without losing any precision during the computation. 

In this paper, we are going to describe two general methods from which one can get the 
integer wavelet transform desired. All of the wavelet transforms from the methods given in this 
paper possess the property of precision preservation (PPP). We draw on the work of several other 
authors who have already contributed to this area [2-3], where some specific examples were 
developed. However, this paper presents a more general method which allows one to see several 
new results as well as those presented and acknowledged prior to this work. 

This paper is organized as follows: Section 2 and 3 give some examples of integer wavelet 
transforms. The examples in section 2 are the starting point for our approach, and the examples in 
Section 3 show the steps and motivation of our general method. Section 4 indicates how one can 
use the lifting technique to create an integer biorthogonal wavelet transform. The Correction 
technique to generate more general integer wavelet transforms is described in Section 5. Section 6 
describes how to process boundaries in order to apply the integer calculation in finite sized images 
or signals. In Section 7, we prove the integer wavelet transforms developed by both the lifting and 
correction method possess the property of precision preservation (PPP). Some example images 
are also shown in this section. The last section, Section 8, provides the conclusion to this paper. 
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2. Basic Integer wavelet transformations 

We provide the following two examples as the starting point for our new method. For the 
sake of convenience, length, and simplicity, we only discuss the algorithm for a one level 
decomposition and reconstruction and only for a one dimensional signal The extension to two 
dimensions is immediate as the rows and columns can be treated into a sequence of one 

c° m j is the original signal where 
the superscript indicates level and the subscript indicates a particular point in the signal. Also, 
{ c i Xlo* m ^ fe* l£o * m * te d^ropostoon parts at the first level. Here 



It. 



^ . * , if N is an even number, 
1 1 if N is an odd number ; 



{ c » mc * its ^ ow frequency (/) part and high frequency {h) part, respectively. For 

multi-levels, we just treat {cj}^ 1 as {cj}*"p and repeat the procedure again. 



E«mpi*i! A (2,2)«wavelet transform by integer calculation. 

This transformation is similar to a variation of the Haar wavelet transform which uses low 
and high pass analysis (decomposition ) filters given as: 



(1) Compute 

(2) Compute 



n 


0 


1 


\ 


1/2 


1/2 


8. 


1/2 


-1/2 



c° ^ , if N is an odd number. 

Here, Int(x) is an arbitrary rounding function which may have different interpretations. 
For example, lnt(x) can be the integer which is nearest to x , or Int(z) may be any integer 

which satisfies x - 1 < Int(jr) S x , etc. It is easy to see that all entries in both {c\ and 



UN is an even number, 



(2.1) 



(2.2) 
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From (2. 1X22), we can easily get the following integer reconstruction algorithm: 



(b) Reconstruction 

( 1) If N is an even number, compute 



or, if n is an odd number, we have 



ck,=ci-lnt(i|o. * «<)...■,*, -i ; (23) 



(2) Compute 



*u.x = c »- InK-^-) , * = 0,- • N t - 2, 



<& m 4+c 9 UM , k = o, -,M,-i. (25) 



Book Since (2.1H2.6) are not linear because of the rounding operation Int(x) , this 

mean, the ^formation order becomes signifies, For instance, if the composition 

waswhedf^to^^ 

appued first to the rows and then to the columns. 

faamplt g; Lazy wavelet transform. 

The Lary wavelet transform does no, do anything. However, this illustrates an important 
concept. The corresponding inverse transform is nothing else but sub-sampling the even 
and odd indexed samples. Decomposition and reconstruction can use same formula as 
follows: 

. < *c«. ft«o,-..,Ar,-t 
d l -cJL» * = <),•..,*/,- l 

m h ^ amPl " 1 2 ^ ^ 800d lMnifoBm for iraa « c but they are simple. 

IlTas1^Tr beaCWeVedfr0n,thCSClWa ^^csteoabove.weco.idTLm 
only as a starting point for our integer, reversible, wavelet transform algorithm 

We ^n™ti°n that there 
which may not be easily seen. If the values of the signal pixels are represented by a finite number 

U» forward transform within the computer itself because of the complementary code property. 
While. ^from the reconstrucUon algorithm, the computer will get back the exact original signal 

^ 5 *T! a COraPlemCntary COdC FOpCny ' Wc «* P-P^y « »W cf Proton 
Preservation (PPp) for these wavelets. 
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It is known that the general values of the high frequency wavelet coefficients are small, and 
all higher levels in the decomposition also provide generally small values in the high frequency 
band. This allows the preservation of precision during the computational stage of the wavelet 
coefficients. Now, the complementary code property, the other aspect of the PPP property is a 
well known characteristic of integer arithmetic as done by the computer. Consider the computation 
of the difference of two integers given as c = b-a and the inverse computation of a m b • c. The 
nature of the computation within the computer can be specified as follows: 

\b-a if -2'*' Sfr-o<2*" 1 -l 
c m = -2l + b-a if b-a>2"~ ] 
2'+b-a if &-fl<-2 9 "' 



and the inverse is 



2«+*-C 



if -l"- x Zb-c m <2 q - x -\ 
if b-c m *V~ x 
if b-c m <-2<-' 

where the m subscript indicates the internal representation, and the range of the integers a, b, c is 
[-2 r \ 2* A - 1], The internal representation of c m when it is outside the range, its appearance is as 
a two's complement number, so the representation may not be the same as the external 
representation of c. However, the same complementary code for the a m will cause the internal 
representation to be identical to the external representation of a. For example, if we let b=2 
(00000010) and a=-/27 (10000001) then c m has the internal binary value of (10000001) when 
q*4. With a value of -I27 for c m , the inverse value for a m will just be o. 

In fact, for Example 2, this property is obviously true. While for Example l, if the range 
of the pixel values is within a finite number of bits, say q, we can only use q bits as the working 
unit, which means the value of transform coefficients will also be within the interval with length 
2*, say [-2 r \ 2*" J - 1]. Due to the nature of computation on a machine, most machines will 
implement (2.1M2.2) automatically as follows (the complementary code property): 



4 = 



0 0 



if •2 rl ScJ 1 -cl t)l <2' H , 



4-^-2'. 



*2«"\ 



if £-4<<-2« 



(2.6) 
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if -2^SInl(i) + ^<2H 



tat 4-) + 4.,-2'. if Tnt(^) + c»^2«-', 

to^+ej^+y, if i^+c'^-r-v 



(2.7) 



While the reconstruction algorithm (2.3) and (2.5) will be implemented by the computer itself 



r-° - 



as 



c]-Int(|-), 
2«+(cJ-Int(£)j, 

(cJ-Int(^-)j-2«, 



if -r-sci-inK^-xr, 

if c]-Irit(4-)<-2^. 



if cj -Int(^.)>2«" 1 . 



if -2'"Srfi + c^ l <2»- 1 . 
4+cL„+2'. if rfj + -2«-\ 
4 + <W2<. if dj+c^'". 



(2.8) 



(2.9) 



It is obvious that (2.8)-(2.9) are just the reverse of (2.6>(2.7). It is also easy to see that if we 
properly take advantage of the bound in the coefficient sire mentioned above, the algorithm can be 
implemented using a minimal amount of storage. 

3. More Examples and Additional Analysis 

In this section we are going to give more examples which will give some motivation for our new 
approach. 

F i i i mnlr A (2,6>wavelet transform by integer calculation [2J. 



This transfo rmation is similar to using fol 



n 


-2 


•1 


0 


1 


2 


3 




0 


0 


1/2 


1/2 


0 


0 


h 


-1/16 


-1/16 


1/2 


An 


1/16 


1/16 



(a) Decomposition 

Decomposition starts with Example 1 at step (1) and (2), and then upgrades the high 
frequency component at step (3): 
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(1) Compute 

(2) Compute 



4... *-o,. 



(3) Compute 



int(~H cj_, , if Ar i$ an even number, 
o 

■ i f A' is an odd number. 



,c'-c' 



4j-lm<SCSL)-4 



1.0 



and then, if Ms even, calculate 



else, calculate 



(b) Reconstruction 



TTw reconstruction algorithm is identical to the decomposition algorithm, except 
running "backwards". 
(1) Compute 



it is now 



4» s Int(i^£LW ( 



4' = im(^Ll£- w ., * sl ,..,„ ( _ 2 , 
and then, if N is even, calculate 

else, calculate 

(2) If /V is an even number, compute 
o i T 
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or, if N is an odd number, wc have 

0 d 1,9 

^««f*-Int(-i-), *B0,- i yv 1 -2. 

(3) Compute 

We see in step (2>-(3) above, that they axe just the same as shown for the reconstruction of 
we (2«2)-wavelet transform (Example 1). 

EttOBbU: A (l,3)-wavclet transform by integer calculation. 

The following nonlinear transform is a variation of the transform which uses biorthogonal 
analysis filters: * 



(a) Decomposition 



n 


•1 


0 


1 


k 


1 


0 


0 


i. 


1/4 


-1/2 


1/4 



This composition starts with the Lazy wavelet at step ( 1 ) and upgrades the high 
frequency component at step (2): 



(DSet 



(2) UN is an even number, calculate 

Otherwisc,if N is an odd number, calculate 

(b) Reconstruction 

(2) If N is an even number, calculate 



,C,s+C 
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Otherwise* if N is an odd number, calculate 

FT H m P Jr * A (5,3)-wavelet transform by integer calculation. 

This transformation is also similar in function to using the biorthogonal analysis filters. It 
is given by 



n 


-2 


.1 


0 


1 


2 


a. 


.1/8 


1/4 


3/4 


1/4 


-1/6 


g. 


1/4 


-1/2 


1/4 


0 


0 



(a) Decomposition 

This decomposition starts with Example 3 at step (1) and upgrade low frequency 
components at step (2): 

0) Set c[°=c£, *»0.-.W,-l; 

If N is an even number, calculate 

Otherwise, if N is an odd number, calculate 

4=Int(^^)-c*,,. * = 0 ( -,M,-l. 

(2) If N is an even number, compute 

ci = ci°-lnt(^), 

«c;=c{'-Int(i^i). fc = l,..,/V,-2, 
Otherwise, if N is an odd number, calculate 
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.ci=ci'-lnt(^i), -2. 

(b) Reconstruction 

(1) Compute 

=c;+Int(-£), 

4-cJ+lM(&ii), i.i, -2. 
4 

Then, if N is even, calculate 

c£-r=<-,+Int( v ' 2 "■-' ). 

4 

J' 

else calculate c j5_, = cJ, iH + Int(-|-) . 

(2) Compute c£ +l = lm(^^) - d\ , * = 0, - . M, - Z 

Then, if N is even, calculate c%_ { = c£_ 2 - d^., . 

The PPP property for Example 1-2 mentioned at the end of the previous section is also 
applicable for these three examples. It is obvious these three transformations are not really linear, 
but they are similar to the one using the corresponding filters given above. Especially, the filters in 
Example 3 and Example 5 belong to, with minor modification, the group of the best biorthogonal 
filters for image compression according to both our experience and the conclusion of [4], 

Also, from the above three examples, we can note that if we begin with integer (linear or 
nonlinear) wavelet transformations and then use some proper upgrading formulas, we can get 
other, much better integer, wavelet transformations for image compression. Now, the key 
problem is: What kind of deductive formulas should be used? We provide an answer to this 
question in the following two sections, Section 4 and Section 5. 

4. Lifting Scheme and Integer Biorthogonal Filtering 

The Lifting scheme, discovered by Sweldens [I], is a new approach for constructing 
biorthogonal wavelets with compact support. However, the most interesting part of this method 
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for us is; it can be used, with minor modification, to create integer biorthogonal wavelet 
transformations. The following is an adaptation of the technique of [1]. 

Definition 1 . The set of filters {A, A, g, g} is a set of biorthogonal filters if the following formula 
is satisfied; 

VcueR m(a>)m f (a)si. 

g v fh((o) h{co+n) 
where m(G)) = \ . . . x , 

and A(0)=2* 4 r to and g(o)=Zg 4 e" to . 

and similarly for m(o), A(G>) and g(fi)) . 

The following lemma is the main result of the lifting scheme [I] reported as corollary 6 in 
that paper. 

Lemma 1. Take an initial set of finite biorthogonal filters {A, A 0 , g °, g ) , then a new set of 
finite biorthogonal filters ( A, A , g, g } can be found as 

Similarly, if we take {A Q , A , g, g 0 } as an initial set of biorthogonal filters, a new set of finite 
biorthogonal filters {A, h , g, g } can be found as 

h(w) = h\(D)+g(Q>)TQGfi 

(4.1b) 

g(<o)=g\o»-h(o>)'&o)- 
Here, s(oi) is a trigonometric polynomial and the corresponding filter s is finite, and so is s(G>). 

Actually* regarding the filters, (4.1) is equivalent to 

' (4.2a) 

— TS - (4 ' 2b) 

Next, we use the lifting scheme with minor modification to create an integer, nonlinear, 
quasi-biorthogonal, wavelet algorithm. Suppose {c J} is a original signal, {c^}and {d|}are 
again its low and high frequency decomposition parts, obtained by using the filters (A, A , g, g J . 
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If we use filters {h , g ) for decomposition (analysis), the corresponding decomposition algorithm 

A 



15 



While the reconstruction algorithm will be 



related to the synthesis filter {A, g) . Here, parameters o t and or, are positive constants with 
a t* a 4 " 2- F OT example, in the situation of regular biorthogonal decomposition and 
reconstruction, a f = a, = ^2 ; and for Example 1 through Example 5 above, a, = 1 and a, = 2 . 

If the set of filters {A, A, g, |} is from {fc, h°, g°, g) by (4.2b), then the decomposition 
can be accomplished as follows: 



Calculate 



Calculate 



The relative reconstruction scheme wil] be: 
1 . Calculate 



/■'•° _i±_V . 



(4.4) 



(4.5) 



2. Calculate 



•"-2yf c «' 0/, "-» i d >£-A 



(4.6) 



Here, equations (4.3) and (4.6) are just the wavelet ( inverse ) transforms using biorthogonal 
filters [h, h 9 , g 9 , g ) . While (4.4) and (4.5) are forward and backward upgrading formulas. 

Similarly, if the set of filters [k, h, g, g) is from the initial set of filters {ft 0 , h , g, g 6 } 
by using (4.2b), the relative decomposition is: 



1. 



Calculate 



4 =a t Ic^..„. 
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2. 



Calculate 



The reconstruction scheme is: 



2. 



1. 



Calculate 



Calculate 



4' = 4+^1^, 



A or, cr. 



For the sake of clarity, we haven't considered the boundary situation, but we will address this 
later. 

Corollary 4. J . Suppose biorthogona! Filters [h, h, g, g) arc from initial filters {h, h\ g\ g) by 
the lifting scheme (4.1a) or (4.2a). If the decomposition and reconstruction by filters 
i*. A 0 . *°, 8 } can be accomplished only by integer calculation, such as Example 2, we also can 
create a corresponding integer wavelet decomposition and reconstruction scheme which is very 
"close" to the original one by using filters {/», h, g, g) . Here the word "close" means that the 
difference of the two decomposition schemes is just some rounding error, and this rounding error 
will be corrected by/the integer reconstruction scheme. 

In fact, if {c| 0 } and {<fj}are integer after (4.3), we can calculate {cj } by 



instead of (4.4). Here Inu» , as described in Section 2. is an arbitrary rounding up function 
which satisfies x - 1 S Int(x) S x + 1 . It is obvious that (4.7) is very "close" to (4 A), and the exact 
reconstruction scheme can easily be obtained from 



and (4.6). Thercwmbeasiinilarresult,ifthesetofbiorthogonal filters [h. h, g, g) is obtained 
from the initial set of filters [h°, h . g, g°) by using (4.2b). 

We can now note, except for the example shown in the Lazy wavelet, (Example 2) most 
standard biorthogonal wavelet transforms cannot be performed directly by integer, even for one of 
the simplest wavelets, the Hoar wavelet. However, if we properly choose the parameters a r and 
or, . and slightly change the transform algorithms, such as Example 1 and Example 3, we can have 
a variation of the original biorthogonal wavelet transforms with respect to the set of filters 




(4.7) 
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{h. g°. g] (or [h Q t h t g, J 0 }). On the other hand, the parameters {j, } should be also 
chosen carefully to guarantee that only addition and shift operations are needed by the algorithm. 

Another observation: if the set of filters {h 9 h , g, g ) is obtained from a set of filters 
{h Q * A. g. g°] by the lifting scheme, and the set {h°, h % g % g 0 } is also obtained from a filter set 

h Q 9 g°, g 0 } , we can repeatedly use Corollary 1 to get a "close" integer wavelet 
transformation. 

5. The Correction Method for Creating Integer Wavelet Transforms 

In this section, we will describe another approach for obtaining integer wavelets by using 
the so called Correction method. The motivation of this method is from the S+P transform, and 
we will now generalize this approach. Actually, the lifting scheme for generating biorthogonal 
wavelets can be considered as a special case of the correction method. From this method we can 
get some even complicated filters with fast decomposition and reconstruction algorithm. 

Suppose that we already have a simple integer wavelet transform, such as Examples 1 
through 3, the decomposition and reconstruction scheme of which can be formulated as follows: 

<f 

Decomposition /r.n\ (5.1) 

Reconstruction c° m = (ft °}{4 L0 J (5.2) 

Here, (5.1) and (5.2) can be the same as (4.3) and (4.6) or other algorithms. 

In general, after the above decomposition, one may not be satisfied with the result. There 
may still be some correlation among the highpass components because of the aliasing from the 
lowpass components, or the lowpass components do not carry enough of the expected information 
from the original signal Hence, we could make an improvement by putting some correction part 
on the highpass components or lowpass components. There are many ways to accomplish this. 
However, for the sake of the integer calculation, we prefer to use following collection method. 
For example, if we want to make a correction for the highpass part, the corresponding formula 
would be: 

d\ sdi-°-Int(<fcl) k=-. 0, 1. 2.-". (5.3) 
Here, dc\ is a correction quantity for d\ 

*i + k « 0, 1 , 2," (5.4) 
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and, {a ( jj^ and |r ; are given parameters which have been chosen for the user's purpose, 

such as reducing the redundancy among bighpass components or some other special requirement 
We are not going to discuss how to choose these parameters, but one can refer to the references [3, 
5, 6) (or clarification of this jwocch, The only thing we need to mention it, for the sake of the 
integer calculation, any entries in both ^ and {r y £ ( should be rational numbers with 

denominators being powers of 2. 

From (5. 1 ), (5.3) and (5.4), it is easy to see the perfect reconstruction algorithm can be 

= ^ + Tnt^, ) k = • f m, m-1, m-2-». (5.5) 

combined with (5.2). 

As mentioned above, the Lifting scheme is a special condition of the Correction method. 
Examples 3 through 5 can also be considered as the examples of this method. We next give an 
example of the Conection method which cannot be included in the group of Lifting scheme, and 
also which does not result in a closed form of compact support for biorthogonal filters* 



K«arwnu /> s+P transform [31, which is similar to using following a nalysis filters 



n 


-2 


-1 


0 


1 


2 


3 


h. 


0 


0 


1/2 


1/2 


0 


0 


X. 


-1/16 


-1/16 


15/32 


-17/32 


7/32 


-1/32 



While, the synthesis filters do not have compact support. However, the S+P transform can 
be implemented as follows: 
(a) Decomposition 

( 1 ) Take the decomposition step of Example 1 , that is, compute 



and 



1.0 



cj =Int<^-)+ c£ u ,. * = O.-.tf, -2. 

1J0 



"M-1 



Int( 1 ) + c HA , if N is an even number. 



, . if N is an odd number. 

(2) Correction Step: Define s 0 5, = 1, 7 *1 and 



and now compute 
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4 

< d\=dr- Int( ^-'* c '- 3c M- 2 ^ ), ik = 1,...)lf 1 .2; 

O 

(b) Reconstruction 

(1) Compute 

^• 0 = ^d+lnt(^). 

(2) If AT is an even number, compute 

4.,-ci-Int(^.). *»ft...,JV,-l 
or, if A/ is an odd number, we have 

cL. s ci-Int(^-), * = -2, 
c° -c 1 

(3) Compute 

4=<+c^, fc = 0.- -,AY,-l. 

6. Boundary Conditioos 

In the previous two sections, we did not show how to get the integer, wavelet transform at 
the boundaries of signals. However, for all of the examples given above, the boundaries have 
been considered. There are two issues dealing with boundary filtering if we use the Lifting scheme 
or the Correction method to generate the integer wavelet transformations. The first is how to 
process the boundaries which occur in the start-up wavelet transformations. The second is how to 
deal with the boundaries in the deductive formula. If the boundaries in the start-up wavelet 
transform have already been established, then those in the upgrading formula are easy to establish. 
In fact, for the Lifting scheme, the boundaries in both steps should be processed in the same way. 
While, for the Correction method, it is easy to see from (5.3)-(5.4) that one has more choices to 
process boundaries in the second step. Therefore, the only thing we need to discuss here is the 
process by which the boundaries in the start up wavelet transformations are established. Assume 
we begin with compact supported biorthogonal wavelets. 
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Suppose the original signal is {c?}_ 0 . For creating integer biorthogonal wavelet 



transformations we can use the following symmetric extension [7]: 

( 1 ) . If current biorthogonal filters have even length, we extend 
the boundaries of the signal as c% = c 4 °_, , k » 1 ,2; • • ; 

(2) . If the filters have odd length, we do the extension as c° t = cf, k « 1 ,2; - . 
Example 1 through 5 use the boundaries give above. In Example 6, the start up wavelet 

transform uses the above boundaries, but in the upgrading step, another boundary filtering is used. 
In addition, for arbitrarily sized images or signals, one can use the same technique which we 
described in the above examples to deal with this condition, 

7. Some Applications 

Before talking about any applications of the integer wavelet transform given above, we first 
prove that a nice property of precision preservation (PPP), which is similar to the one mentioned in 
Section 2, holds for both the Lifting and Correction upgrading technique. This property is very 
important for many applications. 

r.gmmfl 7 1 Suppose that our integer wavelet transform starts with a pair of biorthogonal filters 
with the PPP property discussed in Section 2, that is, (4.3) and (4.6) possess this property. Then, 
the same property will be preserved in the whole algorithm if we adopt the Lifting scheme to be 
the upgrading formula. 

In other words, Lemma 7.1 states if we only use the working units with the same precision as the 
original signal or image to calculate the wavelet transform developed in Section 4, the equations 
(4.8) and (4.6) are still the backward operations of the equations (4.3) and (4.7). 
Proof. Assume that we only use q bits to represent images or signals, say, the range of the pixel 
values is within 2 H -1] . According to the hypothesis of the lemma, the equations (4.3) 

and its inverse (4.6) have the PPP property. Therefore, what we have to verify here is that the 
equation (4.7) and its inverse (4.8) can preserve the same property. We rewrite (4.7) and (4.8) as 
follow: 



CM) 



and its inverse 



j.o ^ A l 



(7.2) 



Here, 
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In fact, tbe quantity b k would have the same value in both (4.7) and (4.8) if we calculate it 
inthesameway. On the other hand, if the working unit for b k is q bits, the machine will give b k 
another value, say b t ( -2-' 55 b k <2^ ). where b k » not equal to b k in the sense of 
mathematics if the value of b, is beyond the interval [-2«-' ( 2 H - 1] . However, £ will be the 
same in both (4.7) and (4.8). Therefore, the machine will automatically implement (7.1) and (7.2) 



and 







if 


.2<H£ C «.o + £ <2 n ( 






c; 0 +fe t -2', 


if 


c\ Q + b t >2«-'. 


(7.1m) 




2< + cl° + b t 


if 


c;-° + i;<-2' 1 . 






c[-b k . 


if 


•2*-' $c[-b k <2"-\ 








if 


c\-b k <-2^, 


(7.2m) 




cl-5-2'. 


if 


c\-b k fcr*. 





It is easy to see that (7.2m) is just the backward operation of (7. 1 m), which provides the evidence 
that the conclusion of this lemma is correct. 

It should be mentioned that the coefficients {c } k } obtained by (4.3) and (7. lm) might not be 
the "real" wavelet coefficients using common sense. However, if we still use the working unit 
with q bits precision at the reconstruction step. (7.2m) and (4.6) will give the exact original signal 
back. On the other hand, the coefficients {cj } still keep the most continuity of the "real" wavelet 
coefficients. Therefore, when we repeat the decomposition step on {cj }, most small coefficients 
in its high frequency part will be almost the same as the "real" coefficients (within some 

rounding error), which allows us to still take advantage of the "real" wavelet transform in image 
compression. 

A similar argument can show the same PPP property will hold for the integer wavelet 
transforms generated by the Correction method in Section 5. 

As we mentioned before, for many applications, tbe lossless image compression is as 
important as lossy compression. The integer wavelet transforms give the opportunity to compress 
without loss. It is also obvious that the integer wavelet algorithms can be used wherever ordinary 
wavelets are used, especially in signal and image compression. However, for most computers, the 
integer wavelet transform is much faster than the ordinary one and it uses much less memory. 
The following are some applications illustrating these types of transforms. 
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Application 1. Lossless image compression 

As mentioned at the beginning of this paper, the integer wavelet transformation 
established by the techniques described in this paper can always be used for lossless image 
compression because of the reversible ability. Especially, we can use the PPP property 
discussed in Lemma 7. 1 . We have used this wavelet lossless technology (WLT) for gray 
scale lossless image compression, and we have tried several images* For most natural 
images, the size of wavelet lossless compressed images is much smaller than corresponding 
GIF images. Figure 1 through 4 give some examples. Figure 1 is a standard image for 
compression, Figure 2 and 4 are X-ray images and Figure 3 is a two-value image but we treat 
it as a 8 bit gray scale image in order to compare with the GIF format In fact, if we convert 
Figure 2 to a binary image, better result can be obtained by the JBIG technique. 




(512*512) 
Figortl. Compression Ratio 

WLT: ).9:l/GlF: 1 05:1 



(512*512) 
Pitur»X Compression Ratio: 

WLT4.5:1/GIF; 2.72:1 



"Visioneer may have come up with on 
against the paper blizzard. . .gets pile: 
way to others throughout your compar 



Figure 3. Compression Ratio: WLT: 20.8/ GIf 17.8 (152x794) 
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Figure 4. Compression Ratio: WLT 3.8M/GIF 1.98:1 (1232*1024) 

Application i Large scale medical image compression 

Usually, 12 bits are used to represent one pixel in medical images. In this situation, the 
values of the pixels vary from 0 to 4095. Such images require careful treatment when a 
transform coding method is used for compression. If we use ordinary biorthogonal wavelets, 
the range of the transform coefficients will expand to [-2 16 , 2 1 6 ] when five levels of transform 
are used Therefore, a longer working unit has to be employed, which consumes significant 
computer resources. However, the integer wavelet technique developed in this paper will solve 
this problem. For example, if we use the transforms given in Example 3, 5 and 6, the values 
of transform coefficients will be limited to the range of [-2 1 *, 2°]. Even if we do not use the 

PPP property for these wavelets, 16 bits for the working unit is sufficient for all computations. 
8. Conclusion 

This paper has shown the processes necessary in order to obtain a non-linear, integer, 
biorthogonal, or non-biorthogonal reversible wavelet transform suitable for signal or image 
processing. We have shown how such a transform can be obtained either using the Lifting 
method, or the Correction method. For example, all interpolation wavelets can be modified to be 
corresponding integer wavelets without loosing any properties of original wavelets. In addition, 
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we have shown under certain conditions, the precision of the transform computation on the 
computer can remain at the same precision of the data, thus reducing the need for additional 
computer memory during the transform computation. These are extremely powerful techniques 
when the target data are large images, or the requirements establish a need for speed. 

Although this paper establishes the structure for the integer transform based upon the 
biorthogonal wavelet or some non-biorthogonal wavelet, we do not imply the examples in this 
paper are necessarily the best wavelets for any particular application. However, we do claim if 
one is going to use such a technique, the ideas suggested in this paper will provide the best 
implementation. 
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Announcing, Lightning Strike 711 Image Compressor (LSIQ version 3.0, a Window! 95 tool that 
compresses still images from 50:1 to 200:1 using INFINlTRON's proprietary wavelet technology. LSIC 
is a versatile, easy to use tool for Web and Graphic designers that can handle a wide variety of digital 
image formats, and it includes filters and convenient web tools* Images can be compressed 3 to 5 times 
more than JPEG, while maintaining similar or better image fidelity. Images can be viewed in 2 to 4 
seconds over the Internet rather than 10 to 20 seconds for images compressed under JPEG* This has 
enormous benefits for reducing bottlenecks on corporate networks and the web, and in addition, requires 
less storage space* 




ligbtnlag Strifct Features 



Compression. Images can be compressed as high as 
200:1 using wavelet technology. 

Comprcs&ioa Control. An EASY mode allows the user 
to compress images with minimal input, requiring only a 
decision between more quality or more compression: An 
ADVANCED mode enables the user to select; I) image 
file sixe, 2) compression ratio, 3) PSNK. or 4) master 
level. Web designers will like the one step process to 
control the size of their image flics, thus insuring the 
speed an image may be viewed oo a browser. 

Non Uniform Compression. Regions of an image can 
be selected for less compression to preserve a higher 
image quality while the rest of the image is compressed 
to the specified compression ratio. In this way important 
parts of a picture maintain crucial details while the over 
all picture file can be made as small as possible 



Post Reconstruction Filters. Filler* arc available to 
enhance tbe reconstructed image. At compression time 
the user can preset a control to have these fillers operate 
automatically during reconstruction. Hie filters include; 
quality improvement, sharpen (edge enhancement), 
smoothing, and brighten. 

Transparencies. Tbe user will have the ability to set 
pixels transparent so that a color in the background 
(already on the page) can be seen through the picture. 
This is useful for creative web site developers. 

Progressive Compression. Ad image can be compressed 
so that when it is viewed it will appear quickly, first wiih 
low resolution, and then progressively building up in 
detail as it is downloaded. This insures the viewer does 
not lose interest while the Image is downloaded 
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Cow, Lightning Strike* 85:1 compression 



MtQtllflllTlll.tlc. 

Founded in 1992* INFINrntON, Inc. is a private 
company ibat specializes in the design and marketing of 
high quality image and video compressira aolutioos for a 
wide array of markets. INFINITRON is based in 
Vancouver, BC with labs in Regina, Saskatchewan and 
Demon Te*as. 
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Applcadra Uimniii Strut Imp Ctapimtoi 

Images on the Web 
Pboto Stock 
Data Warehousing 
Catalogues 
Video Games 

CDs (Encyclopedias, Museums. Science, and Medicine) 
Archives (Art, History, Genealogy) 
Medical Imaging 

rerlirminca 

Encode time typically less than 3 seconds for a 320 X 240 
pixel, 24 bit color image on a 133 MHz. Pentium with 16 MR 
RAM. 

Mlnlmun Rett mmeidftd Sislin 

Windows 95/NT OS 

Pentium 1 00 MHz, 8 MB RAM 

2 MB for program files 

10 MB plus to swap image files 

Jtacnmiinimirapniim 

• Netscape Navigator Plug-in 

• Java Applet 

• ActiveX Control 

• Web Site Image Converter 

• Lightning Strike SDK 

• CiML Banner Generator 




Download a FREE demo version of Lightning Strike Windows Compressor from: 

wwwinfinitronxom 

INFINITRON 3401 East University, #1 04, Denton. TX, 76208 
USA Office Tel: 817.4K4.1165 FAX: 817.484.0588 



10 th Fir 1 199 W. Hastings, Vancouver, BC, V6E 3T5 
Tel: 604.688.9789 FAX: 604.688.9798 



INFINITRON 
Canada Office 

March 1 W7 Crmvrieht INFTNTTRON Research InwmatinnaK Inc. All Riehu Reserved 
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Announcing, Lightning Strike Image Compressor (LSIC) version 3.0, a Windows 95 tool that 
compresses still images from 50:1 to 200:1 using InCnitron's proprietary Wavelet technology. LSIC is a 
versatile, easy to use tool for Web and Graphic designers that can handle a wide variety of digital image 
formats and includes filters, and convenient web tools. Images can be compressed to files 4 times smaller 
than JPEG, while maintaining similar or better image fidelity. Images can be viewed in 2 to 5 seconds 
over the internet rather than 10 to 20 seconds for images compressed under JPEG. This has enormous 
benefits for reducing bottlenecks on corporate networks and the web, and in addition, requires less 
storage space, 




LigMiini Strike Features 



Compression. Images can be compressed to' war 1 
using Wavelet Technology. H p i 

Compression Control. An EASY mode allows the user 
to compress images with minimal input, requiring only a 
decision between more quality or more compression. An 
ADVANCED mode enables the user to select; 1) image 
file size, 2) compression ratio, 3) PSNR, or 4) master 
level. Web designers will like the one step process to 
control the size of their image files, thus insuring the 
speed an image may be viewed on a browser. 

Non Uniform Compression. Regions of an image can 
be selected for less compression to preserve a higher 
image quality while the rest of the image is compressed 
to the specified compression ratio, in this way important 
pans of a picture maintain crucial details while the over 
all picture Hie can be made as small as possible 



Post Reconstruction Filters. Filters are available to 
enhance the reconstructed image. At compression time 
the user can preset a control to bave these filters operate 
automatically during reconstruction. The fillers include; 
quality improvement, sharpen (edge enhancement), 
smoothing, and brighten. 

Transparencies. The user will have the ability to set 
pixels transparent so that a color in the bacVgruund 
(already on the page) can be seen through, the picture. 
This is useful far creative web site developers. 

Progressive Compression. An image can be compressed 
so that when it is viewed it will appear quickly, first with 
low resolution, and then progressively building up in 
detail as it is downloaded. This insures the viewer does 
not loose interest while the image is downloaded. 
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About INFINITRON. Inc. 

Founded in 1992. INFTNTTRON, Inc. is a private 
company that specializes in the design and marketing of 
high quality image and video compression solutions for a 
wide array of markets. INFINITRON is based in 
Vancouver, BC with labs in Regina, Saskatchewan and 
Denton Texas. 
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Itpllctftm Untitling SMki buge Cwnorwsito 

Images on the Web 
Photo Stock, 
Data Warehousing 
Cataloqucs. 
Video Games 

CDs (Encyclopedias, Museums, Science, and Mediuine) 
Archives (Art, History, Geneoloyy) 
Medical Imaging 



Pirlirmiice 

Encode time typically less* than 3 seconds for a 320 x 240 
pixel, 24 bit color image on a 133 MHz. Pentium with 16 MB 
RAM. 

Mlnimun Recinmiidil Syilim 

Windows 95/NT OS 

Pentium 1 00 MHz, g MB RAM 

2 MB for program tiles 

10 MB plus to swap image flics 

Rncllin INHHfTIOR frtducts 

• Netscape Navigator Plug-in 

• Java Applet 

• ActiveX Control 

• Web Site Image Coovener 

• Lightning Strike SDK 

• GMT. Banner Generator 




Download a FREE demo version of Lightning Strike Windows Compressor from: 

www.infinitron.com 

INFINrTRON 3401 East University, #104, Denton, IX, 76208 
USA Office Tel: X17.484.1 165 FAX: 81 7.4 84. 0588 

INFINTTRON 1 0 a Fir 1 1 99 W. Hastings, Vancouver, BC, V6E 3T5 
Canada Offiuc Tel: 604.688.9789 FAX: 604.688.9798 
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The Lightning Strike Compressor i$ a Window* tool that compresses* still images from a wide 
variety ofdigitil image formats using Infinitron's proprietary Wavelet algorithm. Images can 
be compressed to file* 5 times smaller than JPEG, while maintaining simitar or better image 
fidelity, Images can be viewed in 1 or 2 seconds over the internet rather than 10 to 20 
seconds for images compressed under JPEG. This has enormous benefits for transmitting 
over corporate nerworks or the web, and in addition, saves space required for storing all those 
images. 

Lightning Strike is a collection of tools in a user friendly environment Two levels of user 
control are offered, one quick and easy for most applications, the other a master level for the 
advanced user who wishes to control parameters to maximize image quality. 

The compression approach used by Lightning Strike is based upon integer wavelets. This 
technology is acknowledged by leading expert* ax a superior compression technique as 
compared to discrete cosine transform used in JPEG. 



. res";.--;- 

Image Compression Options and Control 

Compression Technique Options 

Both Infinitron's Wavelet Compression and other frequently used compression methods are 
included in the product so users need only have Lightning Strike on their work station to 
perform all image compressions. Images can be compressed to Wavelet, JPEG, PNG, and 
GTF. 

Compression Quality Versus S peed Options 

The user can select one of two encoding processes that trade quality for speed of compression 
and ease of use. With the "Advanced- option selected, the optimum compression parameters 
are set by tho user to give the best possible images for selected compression ratio. With 
"Easy" selected you get the fastest compression without having to know details of parameter 
selection. 

Compression Ratio Control 

The compressed image file size or compression ratio may be specified rather than the quality 
ftctor. This enables a web designer to control the size of their image files or the speed an 
image may be viewed, in a one step process. 

Region of Interest Focusing 

Regions of an image can be selected for less compression to preserve a higher image quality 
while the rest of the image is compressed to the specified compression ratio. In this way 
important parts of a picture maintain crucial details while the over all picture file can be made 
as small as possible. This is also known as Non-Uniform Compression. 
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Spirt and Merge 

Very large images, which could not otherwise be compressed due to their large sire, can be 
split into smaller images and compressed individually. This process has the side advantage of 
using RAM more effectively speeding time for compression. The split images can be 
reassembled using the merge aspect of the feature. 

Post Reconstruction Filters 

Filters arc available to enhance the reconstructed image. At compression time the user can 
preset a control to have these filters operate automatically during reconstruction. The filters 
include; quality improvement, sharpen (edge enhancement), smoothing, and brighten. 

Transparencies 

The user will have the ability to set pixels transparent so that a color in the background 
(already on the page) can be seen through the picture. This is useful for creative web site 
developers. Tms gives the ability to display pictures other than the rectangular shape allotted 
on the web page, i.e. circles, polygons etc.. Also, designers often use this feature for 
shadowing, letters and objects. 

Progressive Decompression 

An image can be compressed so that when it js viewed it will appear quickly, first with low 
resolution, and then progressively building up in detail as it is downloaded. This insures the 
viewer does not loose interest while the image is downloaded. 
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PNG File Structure 

The compressed images tre stored in file format compliant with the PNG standard. In the 
future this file format will replace the GIF format used today. 

Batch Compress 

The user is able to compress many images at once by adding or deleting image file* (or paths) 
to a list box. 

Image Statistics 

The compressor stores image statistics on each compressed image which may be viewed by 
the user. The following information is provided; image dimensions, compression ratio, file 
sizes, MSE, PSNR, maximum pixel difference, compression and decompression times. 



Encoding and Decode Time. 

The typical time to encode or decode a 320X240, 24 bit color image is 1 second on a Pentium 
running at 133 MHz with 16 Meg RAM. 

Minimum Recommended System 

The minimum system requirements for an IBM PC Compatible are: 

Hard Disc Drive 2 Mbytes free for program files, 10 Mb plus to swap 
image files. 

Operating System MS Windows 3.1 (Win32y 95/ NT 
RAM 8 Mbytes 

This software is also available on the Apple MAC, Solaris, and UNIX platforms. 
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Netscape Navigator Plug-In 

Netscape plug-ins are available for (he Mac68k. Mac PPC. Widows 3. 1 and 95/NT. 
Java Applat 

Java Applet* are available for the Mac68k, Mac PPC, Widows 3, 1 and 95/NT. 
ActiveX Control 

The Lightning Strike decompression software is available for such applications as Microsoft's 
Internet Explorer, as an ActiveX control. 

Web Sftt Image Converter 

Thia utility will automate the conversion of web pages from JPEG to the Lightning Strike 
format The utility searches an HTML tile and replicates it replacing any JPEG image tag* 
with Lightning Strike tags and converting the JPEG image Hies to Lightning Strike. The 
utility can follow link tags to recursively convert and replicate an entire web site or sub- 
section of a web site to the Lightning Strike format Thia utility will be available for 
Windows NT and most flavors of the UNIX operating system. 

Lightning Strike Software Developers Kit 

Using the SDK, a developer can integrate the highly efficient Lightning Strike module 
libraries into their own applications. 




Download a FREE demo version of Lightning Strike Windows 
Compressor from: www.infinitron.com 

Infinitron 3401 East University, #104, Denton, TX, 76208 
USA Office Tel: 817.484.1 165 FAX: 817.484.0588 
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I N F I N I T R O N Windows Image 

Compressor V3.0 

Announcing, Lightning Strike™ Image Compressor (LSIC) version 3.0, a Windows 95 too! that 
compresses still images from 20:1 to 200:1 using INFINITRON's proprietary wavelet technology. LSIC 
is a versatile, easy to use tool for Web and Graphic designers that can handle a wide variety of digital 
image formats, and it includes fdters and convenient web tools. Images can be compressed at ratios well 
in excess of present JPEG ratios while maintaining comparable image fidelity. This translates to much 
shorter image down load time on the web. This has enormous benefits for reducing bottlenecks on 
corporate networks and the web, and in addition, requires less storage space. 
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Compression. Uses a proprietary integer wavelet 

Compression Control. An EASY mode allows the user 
to compress images with minimal input, requiring only a 
decision between quality and compression. An 
ADVANCED mode enables the user to select; 1) image 
file size, 2) compression ratio, 3) PSNR, or 4) master 
level for professionals where every parameter can be 
altered. We also provide the highest, wavelet lossless 
compression for users wishing this capability. Web 
designers will like the one step process to control the size 
of their image files allowing control over the delivery 
time of an image over a network. 

Non Uniform Compression. Regions can be selected for 
less compression to preserve image quality while the rest 
of the image is compressed to the specified compression 
ratio.. In this way, important parts of a picture maintain 
crucial details while the over all picture file can be made 



as small as possible 

Post Reconstruction Filters. Filters are available to 
enhance the reconstructed image. At compression time 
the user can preset a control to have these filters operate 
automatically during reconstruction. The filters include; 
visual quality improvement, sharpen (edge 
enhancement), smoothing, and brighten. 

Transparencies. The user will have the ability to set 
pixels transparent, so that a color in the background 
(already on the page) can be seen through the picture. 
This is useful for creative web site developers. 

Progressive Compression. An image can be compressed 
so that when it is viewed it will appear quickly, first with 
low resolution, and then progressively building up in 
detail as it is downloaded. This insures the viewer does 
not lose interest while the image is downloaded. 
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About MFIIITIOI, lie. 

Founded in 1992, INFINITRON, Inc. is a private 
company that specializes in the design and marketing of 
high quality image and video compression solutions for a 
wide array of markets. INFINITRON is based in 
Vancouver, BC with offices in Regina, Saskatchewan and 
Denton Texas. 
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• Images on the Web 

• Photo Stock 

• Data Warehousing 

• Catalogues 

• Video Games 

• CDs (Encyclopedias, Museums, Science, and Medicine) 

• Archives (Art, History, Genealogy) 

• Medical Imaging 

Ptrlirnioct 

Encode time typically less than 2.5 seconds for a 640 X 480 
pixel, 24 bit color image on a 133 MHz. Pentium with 16 MB 
RAM. Decode time is less than .75 seconds. 

MfQlmum licimniaditf Systin 

Windows 95/NT OS 

Pentium 100 MHz, 8 MB RAM 

2 MB for program files 

10 MB plus to swap image files 

Ancillary IMHMmoM Pradocti 

• Netscape Navigator Plug-in 

• Java Applet 

• ActiveX Control 

• Web Site Image Converter 

• Lightning Strike SDK 

• GML Banner Animation/Compression 

• Black and White Image Compression 




Download a FREE demo version of Lightning Strike Windows Compressor from: 

www.infinitron.com 

INFINITRON 3401 East University, #104, Denton, TX, 76208 
USA Office Tel: 817.484.1165 FAX: 817.484.0586 

INFINITRON 10* Fir 1 199 W. Hastings, Vancouver, BC, V6E 3T5 
Canada Office Tel: 604.688.9789 FAX: 604 688.9798 



